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Abstract 

We study the structure of symmetric vortices in a Ginzburg-Landau model based 
on S. C. Zhang's SO(5) theory of high temperature superconductivity and antiferro- 
magnetism. We consider both a full Ginzburg-Landau theory (with Ginzburg-Landau 
scaling parameter k < oo) and a k — ► oo limiting model. In all cases we find that the 
usual superconducting vortices (with normal phase in the central core region) become 
unstable (not energy minimizing) when the chemical potential crosses a threshold level, 
giving rise to a new vortex profile with antiferromagnetic ordering in the core region. 
We show that this phase transition in the cores is due to a bifurcation from a simple 
eigenvalue of the linearized equations. In the limiting large k model we prove that the 
antiferromagnetic core solutions are always nondegenerate local energy minimizers and 
prove an exact multiplicity result for physically relevent solutions. 
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1 Introduction 



In 1986 Bednorz and Miiller announced their discovery of high critical-temperature (To) 
superconductors, and promptly received the 1987 Nobel Prize for their efforts. This dis- 
covery has led to a new flowering of superconductivity theory, since the high temperature 
phenomenon cannot be explained by the accepted models for conventional superconductors. 
In particular, many physicists have come to the conclusion that the microscopic BCS theory 
does not correctly describe the interactions which produce superconductivity at high temper- 
atures. At the present time, there are several competing theories which attempt to explain 
these interactions. One theory is based on the observation that high-Tc compounds also 
exhibit an ordered phase called antiferromagnetism when physical parameters (such as tem- 
perature, chemical potential or "doping", and magnetic field) are varied. Antiferromagetism 
(abbreviated AF) is an insulating phase of matter in which electron spins orient themselves 
in the direction opposite to their nearest neighbors. The coexistence of these two phases 
(AF and SC) in the phase diagram of the high-Tc compounds has led to the speculation 
that high temperature superconductivity and antiferromagnetism could be explained by the 
same type of interaction. 



Following in this direction, Shou-Cheng Zhang ||Z 97]| proposed a quantum statistical 
mechanics model which incorporates AF and high temperature superconductivity (SC). The 
model is based on a broken 5*0(5) symmetry tying the complex order parameter of supercon- 
ductivity to the Neel vector which describes antiferromagnetism. The interactions between 
the SC and AF order parameters in this model should have some effect on the familiar con- 
structions from conventional superconductivity theory. In a recent paper Arovas, Berlinsky, 
Kallin, & Zhang [ |ABKZ 9"7| introduced a phenomenological Ginzburg-Landau model based 
on the 50 (5) theory, and studied isolated vortex solutions in the plane. Recall that in a 
conventional superconductor the magnetic field is expelled from the superconducting bulk, 
and only penetrates in thin tubes (the vortices) where superconductivity is supressed. Hence, 
in the conventional theory the magnetic field is constrained to a small core of normal (non- 
SC) phase. Using a simplified model Arovas et al predicted a new kind of vortex structure 
in the 5*0(5) model: vortices with antiferromagnetic cores, which should be observed for 
small values of the chemical potential. They also predicted that (as the chemical potential 
is gradually decreased) the transition from normal core to AF core vortices occurs in a dis- 
continuous fashion. In other words, AF cores should be produced via a first order phase 
transition. 



In this paper we rigorously analyse vortex cores in the full 50(5) Ginzburg-Landau model 
and in an "extreme type II" limiting model (also called "high kappa model" ) to understand 
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the nature of the transition between normal core and AF core solutions. For both models 
we show that the vortex solutions with normal cores become unstable (within the class of 
radial functions- see ( p. . 1|) below,) and vortices with AF cores are produced by bifurcation 
from the normal core solutions. In the extreme type II model we prove that the transition 
is continuous (ie, second order), contrary to the prediction of |[ABKZ 97|| (see Figure 1.) 
Furthermore, we show that for each value of the chemical potential there exists a unique 
stable vortex profile (see Theorem |4.5| .) 



The full SO(5) Ginzburg-Landau free energy is written in terms of the SC order parame- 
ter ip G C and the AF order parameter (Neel vector) m = (m 1 ,m 2 ,m 3 ). In non-dimensional 
form, the free energy is: 




m\ 2 ) 2 + gK 2 \m\ 2 + |(-V - A)i)\ 2 + |Vm| 2 + |V x A\ 2 1 dx. 



(We refer to the paper by Alama, Berlinsky, Bronsard & Giorgi ||ABBG 98| where the free 
energy is written in dimensional form.) In these variables, the penetration depth A = 1, 
and the Ginzburg-Landau parameter k is the reciprocal of the correlation length £. The 
parameter g measures the strength of doping (chemical potential) of the material. It is this 
term which breaks the SO(5) symmetry of the potential term. We take g > 0: with this 
assumption superconductivity is preferred in the bulk of the sample. 

M 2 we seek critical points of T of 



To study isolated vortex solutions in the plane f2 
the form 

Jd0 A a/„\ (~ V X 



ip = f( r y 



A = S(r) 



m = mirirrio 



where mo a fixed unit vector, and d G Z \ {0} represents the degree of the vortex. As for 
conventional SC vortices, we expect that only the solutions with d = ±1 will be energy 
minimizers (see Gustafson [ |Gu 98[ j , Ovchinnikov & Sigal JUS 97 
this ansatz solve the system of equations 



Critical points of T with 



(GL) K , : 



-f" - If 



i fl , (d-Sf 



-S" + }S' = (d-S)f 2 , 
—m" — -w! + K 2 gm = k 2 (1 — f 2 — m 2 )m, 
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with f{r) > 0, f(r), S(r) — > as r — > 0, and f(r) — > 1; 5(r) — > d as r — >■ oo; and m'(0) = 0, 
m(r) — > as r — > oo. 

In addition, we study the following "extreme Type II" model, 



(GLjoo^ 



-f" ~ If + £/ = (1 - f 2 ~ m 2 )f, 



-m 



-m' + gm = (1 — f 2 — m 2 )m. 



The system (GL)^ ^ is obtained in the limit k — > oo after rescaling solutions to (GL) Kj5 by 
the correlation length £ = 1/k. For high Tc superconductors k is very large, and hence 
the vortex cores are very narrow compared to the penetration depth, which measures the 
length scale for magnetic fields. By rescaling we capture the structure of the vortex cores 
and decouple the magnetic field, which lives on a much larger length scale. Indeed, the 
calculations which led Arovas et al JABKZ W? \ to predict AF vortex cores are mostly based 
on (GL)oo iS and its associated free energy functional. 

We observe that when the AF order parameter m = the two systems (GL) re>g and 
(GL)oo i9 reduce to the familiar Ginzburg-Landau vortex equations, well studied in the math- 
ematical literature (see Plohr ]F~gQj , Berger and Chen | |BG 89|| , Chen, Elliot, & Qi | |CFQ 5g , 
Brezis, Merle, & Riviere [ BMK 94 1, Ovchinnikov & Sigal [ L)S 97 |, for example.) We call these 
the normal core solutions. In a previous paper |[ABG 99| we have proven that when k 2 > 2d 2 
there is a unique normal core solution, which is a non-degenerate minimizer of the appro- 
priate free energy functional. This characterization will be essential for our analysis of the 
normal-to-AF core transition. 



We now discuss our results. We define a reduced energy functional defined for functions 
satisfying the symmetric vortex ansatz ( [Lip , as well as appropriate function spaces in which 
that functional is smooth. We find that for every k (including the extreme type II model) 
there exists g* > such that the conventional normal core vortex solutions of (GL) K)£ , (and 
(GL^oog) are strict local minimizers of the reduced energy for g > g*, but are not local 
minimizers when < g < g*. In particular, energy minimizers must have AF order in the 
vortex core for < g < g*. When k 2 > 2d 2 we show that the AF core solutions bifurcate 
from the normal core solution at a simple eigenvalue of the linearized system (GL) K;5 (or 
(GL)^^.) The bifurcating solutions remain bounded for g > and lose compactness as 
g — > 0+ with / — > and m — > 1. 

For the limiting problem (GL)^ we obtain a complete picture of the phase transition 
to AF cores. This is because all AF core vortex solutions are non-degenerate minima of the 
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reduced energy. (See Theorem |3.1| .) Stable (locally minimizing) solutions with m(r) > 
bifurcate from m = at g = g*^ to values g < g^. Moreover, for each g < g^ there exists 
exactly one solution with m(r) > 0. 

In the language of physics, our results indicate a second order (or continuous) phase 
transition between normal and AF vortex cores in (GL)oo i9 . This information concerning the 
nature of the transition was not derived in the paper by Arovas et al |[ABKZ 97| , and hence 
the result is new to the physics literature as well. For (GL) KjS Alama, Berlinsky, Bronsard, & 
Giorgi [|ABBG 98|| present numerical simulations (based on gradient flow for a finite elements 
approximation of the free energy) which suggest that the transition is also second order for 
k < oo. (See Figure 1.) However we were not able to extend the arguments used in studying 
the bifurcation curves of (GL)^ ^ to the more complicated system (GL) K g . See Remark |4.3| 
for further discussion. 

Here is an outline of the content of the paper. In the second section we introduce the 
reduced energy and function spaces, we treat briefly the questions of existence, regularity, and 
decay of solutions, and we present properties of physically relevant ( "admissible" ) solutions. 
We also prove the monotonicity of the solution profiles (/, S, m) under the hypothesis that 
the solution is a local reduced energy minimizer. This result (Theorem |2l|) is done in the 
spirit of the weak maximum principle (see Theorem 8.1 of ||GT 83|| .) 

Section 3 contains the proof that all solutions of (GL)^^ with m > represent non- 
degenerate local minima of the reduced energy. This result is the key to understanding the 
bifucation diagram for (GL)^^. The bifurcation analysis itself occupies Section 4. 

The last two sections contain the a priori estimates used in rigorously passing to the limit 
k — ► oo and in studying the global behavior of bifurcating continua. In both cases, we require 
estimates on solutions which are energy- independent. For the limit k — > oo this is because the 
reduced energy of minimizers behaves like log k, and in studying global bifurcation we require 
estimates valid for any physically relevant solution (whether it is energy minimizing or not.) 



The starting point for these estimates is a Pohozaev type identity (see Proposition |5. 4j .) The 
proof of convergence to (GL)^^ as n — > oo is presented in Section 5; other a priori estimates 
are derived in Section 6. 

We wish to thank our colleague John Berlinsky for introducing us to the 5*0(5) model, 
and for his great patience in explaining physics to we mathematicians. We are also obliged 
to the Brockhouse Institute for Materials Research for supporting a workshop which brought 
together physicists and mathematicians to discuss issues in superconductivity. 
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k=20 (dotted), 40 (dot-dash), 120 (solid) 
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Figure 1: A numerical bifurcation curve, m(0) vs. g for values k = 20,40, 120 and d = 1, 
indicates a second-order transition to AF cores in model (GL) K g . For the high-kappa model 
(GL^oo^ we prove that the above image correctly depicts the solution set (see Theorem [15].) 
Numerical simulations indicate that the bifurcation occurs at ~ 0.2545 ||ABBG 98| . 



S. Alama, L. Bronsard, T. Giorgi 



6 



2 Solutions of the Ginzburg— Landau system 
2.1 Preliminaries 

Here and in the rest of the paper, we fix the value of d £ Z \ {0}. In this section k £ M 
is fixed. Note that without loss of generality we may take d > 0, since the free energy 
and the corresponding Euler-Lagrange equations are invariant under the transformation 
(V>, A,m) -> (i>, -A,m). 

Following our previous work ||ABG 99]j on symmetric vortices, we define a function space 
for which the free energy will be a smooth functional. First we fix some notation: we denote 
by LP, H the Lebesgue and Sobolev spaces (respectively) of radially symmetric functions in 
J? 2 , that is, 



rod 

L v r = {u(r) : / \u(r) \ p r dr < oo}, (p<oo), 
Jo 



H:=Hl = Mr): \(u'(r)) 2 + (u(r)) 2 ] r dr < oo}, 

jo L J 



and analogously for Lf. We also denote / u{r)rdr = f£°u(r)rdr. 
Define the Hilbert space 



f u 2 

X = {u £ H : —rdr< oo}, 
J r 2 



with norm 



\u\\x 



\ 



(u'(r)) 2 + u 2 + — 



r dr. 



The following density and imbedding properties for the space X are proven in [|ABG 99|| : 



Lemma 2.1 i. X is compactly embedded in L v r for each p £ (2, oo). 
it. X is compactly embedded in L 2 rloc . 
Hi. For every u £ X, 



Woo < 



u'? + 



u 



r dr. 
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In particular, X embeds continuously into L£°. 
iv. C£°((0, oo)) is dense in X. 

We note that the compactness of the embedding of H into L p rloc (1 < p < oo) is just 
the classical Rellich-Kondrachov Theorem, and the compact embedding of H into LP r for 
2 < p < oo is due to Strauss [Bt 77 . 



2.2 Energy 

We now define our energy functionals, using the space X defined above. To keep the appro- 
priate boundary condition at infinity we fix any function i] G C°°([0, oo)) with rj(r) = for 
< r < 1, r](r) = 1 for all r > 2, and < r\ < 1. Then set /o = r/, So = drj, and seek 
solutions (/, S,m) of (GL) Kj9 with f = f + u, S = S + rv, u,v & X, m & H. (Later we 
will see that this choice poses no restriction on solutions which are physically relevant.) We 
denote by Y = X x X x H, and by Y the affine space 

Y = {(f,S,m) : / = /„ + «, S = S + rv, u,v e X, m G H} = Y + (/„, S , 0). 

For (f,S,m) G K we define 
(2.1) £ K>g (f,S,m) = 



U'Y + 



5_' 

r 



1 2 



/ / \ 2 i 2 2 

(mj + k gm 



+ {d r f 2/2 + y(l-/ 2 -^ 2 ) 2 | 



r 



and the functional Jfc iS : Y — > iR by 

4, g (tt, u, m) = £ Kt g(f + u,S + rv, m) - £ K , g (fo, S , 0). 

Throughout the paper we will take advantage of these two representations of our spaces and 
energies, and use the formulation which is more convenient at the given moment. 

Defining an energy functional for the limiting problem (GL)^^ is trickier, since the naive 
choice for the energy (namely (|2.1|) with 5 = and k — 1) would be infinite for all / 
satisfying the desired boundary condition at r = oo. Our solution is to subtract off the 
offending term from the energy density. Let be the (unique) positive solution to the high 
kappa vortex equation, 

1 ~/ d 2 ~ ~ 2 ~ 

~foo foo ~l of <x> ~ (1 — foo)foo 
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with /oo(0) = 0, foo(r) — > 1 as r — > oo. The uniqueness of f^ was established by Chen, 
Elliot & Qi [CEQ 94| |. The estimates in [ CEQ 94| ensure that is smooth, /oo(r) ~ r d 
near r = 0, and (1 — foo) G 7F 

We define the appropriate spaces for the free energy Soo, g based on f^: let Z Q = X x H 
and 

Z = {(/,m): f = L + u, ueX, m G H} = Z + 0). 
Then the energy for the high kappa model is: 

(2.2) W/>™) = 



(/' ) 2 + K) 2 + am 2 + ^ [/ 2 - £] + ^(1 - / 2 - ™ 2 ) 2 i> ' <!>■ 



1 

2< 



If we write / = foo + u, we reduce to the equivalent functional 

(2.3) Ioo,g{u, m) = £oo,g(/oo + % m) - £oo,g(/oo, 0) 



+^(1 - (/oo + n) 2 - m 2 ) 2 - 1(1 - /^) 2 + 2(1 - pJUu \ r dr. 



By a direct expansion of the energy in powers of u,v,m we see that Ik, g : Y ^ 1R and 
Ioo,g '■ Z — > iR are smooth (C°°) functionals. 

When (7 > is fixed, we obtain solutions of (GL) K>g and (GL)^^ as global minimizers for 
S^g and £oo jff (in the appropriate spaces, Y and Z): 

Theorem 2.2 For every fixed g > 0, k G M, d G Z — 0, £/ie functional J Kj5 admits a 
minimizer (u,v,m) G X x X x if. Moreover, (/, S 1 , to) = (/o + u, So + rt>, to) is a smooth 
solution of the system (GL) K ^ g . 



Theorem 2.3 For even/ fixed g > and d G Z — 0, the functional I^g admits a minimizer 
(u,m) G X xH. Moreover, (f,m) = (foo+u,m) is a smooth solution of the system (GL)oo, g - 



The proofs of Theorems |2.2| and Theorem |2.3j are straightforward but technical, and are 
deferred to Section 6. 
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2.3 Admissible solutions 

As in [|ABG 99j], we define a natural class of solutions to the system (GL) KiS 



Definition 2.4 We call (/*, S*,m*) an admissible solution to (GL) K ^ g if: 
i. ( GL) K ^ g holds for all r G (0, oo); 

Hi. f*(r) > and m*(r) > for all r > 0; 
iv. 5,(0) = and r<(0) = 0. 

A solution (/*,m*) of (GL)^^ is called admissible if the above conditions hold, where we 
replace k by oo and disregard S*. 

A solution to (GL) K g or (GLj^ g with m* = is called a normal core solution. 

The admissible solutions are those which are physically relevant in the context of the vortex 
core problem described in the introduction. We note that the normal core solutions are 



unique for k 2 > 2d 2 : see |[ABG 99|| for the case 2d < k < oo and [|CKQ 94|| for k = oo. 



We now present some properties of admissible solutions. In the following, we will assume 
that k G M[j{oo}, with the understanding that S** = when k = oo. 

Proposition 2.5 Let (/*, S*, m*) be any admissible solution of (GL) K g . Then: 

i. For all r G (0, oo) it holds < /*(r) < 1, < m#(r) < 1, f 2 (r) + m 2 {r) < I, and, if 
k 7^ oo, < 5*(r) < d. 

ii. Either m*{r) > for all r G [0, oo), or m* vanishes identically. 

Hi. /*(r) — > 1, m*(r) — > 0, and, if n ^ oo, S^r) — >• d as r — > oo. Moreover, there exist 
constants a, Cq > stzc/i t/iat /or k ^ oo 



< 1 - /„(r) < C e- CTr , < d - S*{r) < C Q e' ar < m,(r) < C e 
and for k = oo 



< 1 - f.(r) + M2 + f + 0(r- 6 ), < m*(r) < C e 



■err 



2r 2 8r 4 

/or all r > 0. 
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iv. /*(r) ~ r d , S*(r) ~ r 2 for r ~ 0. 

v. If k oo, Sl(r) > /or all r > 0. 



Proof: The proof is very similar to that of Proposition 2.3 of | |ABG 99] ], so we provide only 
a sketch. From the fmiteness of the free energy we immediately conclude that m* G H, and 
hence m* G for any p G [2, oo], and ra*(r) — * as r — > oo. Since /* > 0, fmiteness of 
energy again implies 1—/* G L 2 . (see ( |6.10|) for details,) and therefore the bound < f*(r) < 1 



follows exactly as in Proposition 2.3 of [|ABG 99(1 . When k < oo, the bound < S*(r) < d 
and the proof that S'^(r) > are also unchanged from ||ABG 99fl . To show z = f%+m% < 1 we 
use the equation satified by z: this argument is already presented in ||ABBG 98(] . Statement 
(ii) is a simple consequence of the strong maximum principle. 

The exponential decay in (iii) for is consequence of Proposition 7.4 in Jaffe Sz Taubes 
(JT 80|| , and so are the ones for /* and S 1 * if k ^ oo. If k = oo, the polynomial decay of /* 
can be proven as in Lemma 3.3 in ||GEQ 94(1 , since m*{r) < for any r > R with C(R) a 
big enough constant. 



The behavior at zero given in (iv) can be proven as in [P 80 







We now connect admissible solutions to our space X. 



Proposition 2.6 Let (/o, So, mo), (/i, Si, mi) be admissible solutions to ( GL) K ^ 9 . Then (f\ — 
fo) e X, [(5*1 - S )/r) e X, and m b m G H. 

Proof: As already remarked, condition (ii) of the definition of admissible solutions implies 
mi, mo G H, and mi, mo G LP for any p G [2,oo). Then, the rest of the proposition for 
k 7^ oo is proven as in Proposition 2.4 of |[ABG 99(1 • When k = oo, we note that (1 — fj) G H 
for £ = 1,2, and that by (iv) of Proposition 2.5 we have (f\ — ^) 2 — cr 2d for r ~ and again 
by fmiteness of energy we conclude our statement. 
<> 



Remark 2.7 In light of Proposition |2.6| we observe that the choice of fo, So in the definition 
of the space Y may be replaced by any fixed admissible solution of the equations (GL) KjS . 
It will be convenient to choose instead the "basepoint" (f K , S K , 0) to be a "normal core" 
solution to (GL) KjS . In other words, an equivalent definition of the space Y is: 



(2.4) Y = {(f,S,m): f = f K + u, S = S K + rv, u,v G X, m G H} 

We recall that the normal core solutions are uniquely determined for k 2 > 2d 2 . 
k 2 < 2d 2 we fix any one. 



When 
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Remark 2.8 Proposition |2.(j| also implies that the admissible solutions are exactly those 
which arise from minimization problems for S Kt9 and £oo,g m the space Y. In particular, as 
an immediate corollary we obtain the following statement: 

(/*, 5*, m*) is an admissible solution to (GL) K 9 if and only if /* > 0, > 0, (f*,S*, m*) G Y 
and £' Kt9 (f*, S*, m*)[u, v, w] — for all u, v G X, and w <E H. 

An analogous statement holds for the problem (GL)^^. 

With this choice of representation for our spaces Y, Z, we now look at the second variation 
of energy with respect to the variables (u,v,w) G X x X x H. We define 

£",g(f*,S*,m*)[u,v,w] 
(2.5) 



(2.6) S^ g (Um*)[u,w] 



Note that if we write /* = f K + w*, £* = S K + rv *, then 

£", g (f*, S*,m*)[u,v, w] = D 2 I Ktg {u in v, n m je )[u,v, iw], 
the usual second Frechet derivative. 

For admissible solutions which are stable, in the sense that the second variation of energy 
about the solution is a non-negative quadratic form, we have monotonicity of the profiles 
f{r), m{r). 



dt 2 



+ tu, S* + trv, m* + tw) 



t=o 



iu'f + (w 1 ) 2 + {d f )2 



2 , 2 2 i / i\2 
U + K QW + [V J 



v 2 Ad — S*) , 



2„.2 



-k 2 (1 - f 2 -m 2 )(u 2 + w 2 ) + 2K 2 (f*u + m*w) 2 } r dr. 



dt 2 



£oo,g{f* + tU,m* + tw) 



t=0 



d 2 

(it ) +{w) + ^it + gw 



;i - /* - ml)(u 2 + w 2 ) + 2(/,u + m,w) 2 } r dr. 
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Theorem 2.9 Suppose (f,S,m) is an admissible solution of (GL) Kt9 , and £" (f, S,m) > 
as a quadratic form acting on Y . Then f'(r) > and (if it is not identically zero) m'{r) < 
for all r > 0. 

For the problem (GL)^ the same theorem holds, with exactly the same proof. We will see 
later that all admissible solutions of (GL)oo iS with m{r) > are stable (in the above sense), 
and hence we will obtain the stronger result announced in Corollary |3.2| . 

Proof: Let u(r) = f'(r), w(r) = mf(r). Then, differentiating the first and third equations 
of (GL)^, 



_ I# + ( d - k 2 (1 - 3f 2 - m 2 )u + 2K 2 mfw 



1~ n d-S P 
u + 2 / 



S' d-S 
— + — s— 



-w" w + gn 2 w — k 2 (1 — f 2 — 3m 2 )w + 2k 2 fmu = — -w. 



Suppose there exist intervals (a, b), (c,d) such that 

u{r) < r G (a, b), u(a) = = u(b)] or 
w(r) > r G (c,d), w(c) = = w(d). 



Note that by the properties (i), (iii) and (iv) of admissible solutions in Proposition p^5|, a ^ 0, 
b, d < +oo. Let 



u(r) 



u(r), if r G (a, b), 
0, otherwise, 



w(r) 



w(r), if r G (c, d), 
0, otherwise, 



Then u < 0, w > 0, and an integration by parts shows that 



(u) r dr 



b l 

u -(ru'Y r dr. 



a r 



and similarly for w. If we now use (u, 0, w) as a test function in the second variation of 
energy and recall from Proposition |2l| that S(r) < d, S'(r) > for all r > 0, we obtain 



< £" (f,S,m)[u,0,w] 
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1 9 

— + 2- 



d-5 



./ 



"5' 




1 2 


— 






r 




u - -^w 



r dr < 0, 



unless u, w = 0. Consequently, u = f > and w = m' < 0. Strict inequality follows from 
the Strong Maximum Principle, since u, w satisfy equations of the form 

— A r it + Ci(r)u > —2n 2 mfw > 0, 
— A r w + C2(r)w = —2K 2 mfu < 0. 







3 Nondegeneracy of solutions of (GL) 



oo,g 



Theorem 3.1 For any admissible solution (/*,m*) of (GL)^^ with m* > there exists a 
constant o~* > such that 



£Lq(f*, m *) [u,w] > <r*{\\u\\% + \\w\\%), 



for all u G X , w G if. 



Corollary 3.2 For any admissible solution (/*, m#) o/ (GL)^^, /*(r) > /or a// r > 0. If 
m* is not identically zero, then m*(r) < for all r > 0. 



The corollary follows from Theorem |3.1| and the argument of Theorem |2 . 9| when m* > 0. Note 
that when m* = the system (GL)^^ reduces to the single equation studied in [ |CEQ 94j 
and the strict monotonicity of /* is part of their result. Also, in the case that m* = the 
Theorem reduces to S'^ „(/*) [u] > 0"*||«||x- 



The key step in proving Theorem |3J] is the following identity: 



Theorem 3.3 For any admissible solution (/*, m#) of ( GL)^^ with m* > 0, and any u G X , 
w G H, 



u 

7* 



+ 171* 



W 



+ 4(/*w + m^w) 2 > rdr 
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Proof of Theorem |3.3| : First we prove the identity for u G C^°((0, oo)) and w G 
C£°([0, oo)). First, note that using /* > and m* > 0, we have 



(3.2) 



fie- 



u 

7* 



J\2 



[U 



t\2 



U 



f* f* 



2 ' 



with a similar identity holding for m*, w. Hence, 



= Z'ooAf*,™*) 



u 2 w 2 

f* TTl* 



{u'f + {w'f + - 2 u 2 + gw 2 - (1 - fl - m 2 )(u 2 + w 2 ^ 



f * 



u 

J*. 



w 
L\m* 



r dr. 



Substituting this in the formula for S'^ (/*, m*)[it, w] we obtain 



u 

J* 



w 



i 2 



2[/*it + m*u>] 2 > r dr. 



To obtain the result for any (it, w) 6lx H, let it n be a sequence of Cq°((0, oo)) functions 
converging to u in X, and iu n a sequence in Cq°([0, oo)) converging to id in if. By continuity 
of S'^gif*, S*), the limit passes in the second variation of £oo, g - For the right hand side we 
expand, 



>\ 2 



u 

J* 



r dr 



■jr + U. 



(i/) 2 — + { 4r I u 2 \ r dr, 



and note that 



(ff £c(1 4> 
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since /* ~ r d for r ~ 0. Hence each term is controlled by the X-norm and can be passed to 
the limit. A similar argument may be applied for the second term in the right-hand side of 
(|3.1|) . The quotient is expanded as in (§) above, with m^.w replacing Then we claim 

that m'(r)/m(r) is uniformly bounded for r G [0, oo). Indeed, by the basic gradient bound 
for solutions of the Poisson equation (see section 3.4 of [|GT 83|| ) we have for any r > 1, 



|t77. / (t"o) I ^ 2 sup m(r) H — sup |k 2 (1 — g — f 2 — m 2 )m\ < C\ sup m{r) 



\r-r \<l 



\r-r \<l 



\r-ro\<l 



Applying the Harnack inequality (Corollary 9.25 of [pT 83||) we then obtain: 



m'(ro) 



m(r J 



, r SUPl.-^lMO . r SUP|r-r |<l m V) , 

m(r ) mf| r _ ro |<i m(r) 



for all tq > 1. Therefore m'/m is uniformly bounded, and we may pass to the H\ limit in 
the second term in (|3.1|) . The last term is clearly continuous in the L^-norm in both u and 
w. In conclusion, we may pass to the limit u n — > u, w n — > w and obtain ( |3.1|) for u G X, 
w G H. 





Proof of Theorem 3.1: Define 



a* = inf{£^ )g (/*,m*)[«,w] : u G X, w G H, \\u\\ 2 x + \\w\\ 2 H = 1}. 

We must show that a* > 0. 

By Theorem |3.3|, a* > 0. To obtain a contradiction, assume instead that a* = 0. We claim 



that in this case the infimum is attained at a nontrivial (-u*,u>*), with m*)[it*, w*\ 

cr* = 0. But this contradicts Theorem B.3L and hence cr* > 0. 



We now claim that the infimum cr* = is attained in Z . Take any minimizing sequence: 



(u n ,w n ) G X x H with ||M n ||x 



w 



n\\H 



1 and 



0. 



By the Sobolev embedding, there exists a subsequence (still denoted by u n , w n ) and -u* G X, 
io* G if so that M n -> m„, ?«„ -+ w„ weakly in X, H (respectively), and strongly in Lf oc . 
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First, we claim that (u*,w*) ^ (0,0). Indeed, if both u^^w* vanish identically then by 
weak convergence (u n ,w n ) — (u*,w*) = (0,0) and the compact embeddings, 



/ («) 2 + ^ u n + 2fi2 /!«' + K) 2 + 9K 2 w 2 n ) rdr 



£'LAf*i m *)[UmW n ] 



+ 



2/i y-2 2\/ 2 , 2\ o 2 2 2 a 2c 

— > 0. 



r dr 



In particular, (u n ,w n ) — > (0,0) in the norm on I x ff, which contradicts the fact that 
||w n ||x + ll^nlln = 1- Thus the claim holds, and («*,«;*) 7^ (0,0). 

Next, we use lower semicontinuity in the norm and Lf oc convergence to pass to the limit, 



(3.3) 



£'<LM*' m *){ u *' w *] - lin^inf ^,«,(/*> m *)[ M n,w n ] = 0. 



This contradicts Theorem |3l| since „(/*, m*) [it*, to*] > 0. (Note that -u/ /* is non-constant 
since «£l but /* ^ X.) We conclude that a* > 0, as desired. 



We note that the same result holds when m* = 0. Hence following the method of 
[|ABG 99 | 1 , we obtain another proof of uniqueness for the solution to the high kappa equation 
for /, studied in [|CEQ 94 1. 



4 Bifurcation from the normal cores 

In this section we show that (when k 2 > 2d 2 ) AF core solutions are nucleated by means of 
a bifurcation from the normal core solution family at a simple eigenvalue of the linearized 
equations. We will also require a priori estimates (whose proof we will present in Section 6) 
to obtain global information about the solutions set for all k 2 > 2d 2 , and the stronger result 



of Theorem |3.1| to fully categorize solutions in the extreme type-II model (GL)oo )S . We 
present the detailed argument for the problem (GL) K g . The functional analytic framework 
is entirely similar for the problem (GL)^ ^, and so we omit it and concentrate instead on the 
more precise global charaterization of solutions which we prove for (GL)oo )S . 
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4.1 Local bifurcation at g* K 

We define a map T : Y x M — > F * by 

((«, w),F(f*, S*, m*, £?)>y ,y * = m*)[u, w], 

(■u, f , iu) G Yo, (/*, 5*, 777.*) G Y. Its linearization is the operator F'if*, S*, m*, g) G £(Yo, ^o*) 
defined by 

£' K , g (f* + 'S'* + rtip,m* + *£)[ U > U >H- 

t=o 

We remark that the explicit expansion of the energy I Kg in terms of w* = /* — f K , t>* = 
(5* — S K )/r, u>* ensures that T is a C 2 map in all arguments «*, t>*, w*, g. 

By the natural identification Y ~ F * of a Hilbert Space with its dual, we may also 
represent T' by C g G £(F , Y ) as 

((m,w,w),£ 9 [v?,^,^)y = ((u, «,«;), F'(f*, S*,m*, g)[ip,ijj, £]) Yo ,y * ■ 
If « : Z* — > Z is the isomorphism, then C g = i o 5*, m*, o). 

Lemma 4.1 For a// g > 0, C g is a Fredholm operator of index zero. 

Proof: Define an equivalent inner product on Y , 



{(u,v,w),F'(f*,S*,m*,g)[ip,'ilj,t])Y ,Y * = -r 



((u,U,lu), (^,-0,0)^0 = J W^' + 2k 2 U^+ -^—wp + v'^' 



+vip H — -vip + w'C,' + gK 2 w{\ r dr. 



Then we write 



where if is defined by 

((u,v,w),K[cp,tp,^]) Yo = J\^K 2 (fl-l)u(p + 2K 2 f*m*(u£ + w(p) + 2K 2 mlw£ > 
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k 2 (1 - ft - mi)(u<p + <) + (ft - 1)W 



-2- — —f^iiip + vip) 



r dr. 



Recalling the decay properties of /*, £*, m* and the embedding properties of H, X we observe 
that K is compact, and hence C g = Idy + K is Fredholm with index zero. 





As a direct consequence of Lemma |4.1| , 

dimker(jF') = dimker(£ 9 ) = codimRan(£ s ) = codimRan (JF'). 

Now we may apply the standard bifurcation theory of Crandall & Rabinowitz [|CR 71|| at 
an eigenvalue g* of J-'(f K , S K , 0, g*). Indeed, note that when m* = the linearization of T 
decouples into two components, 

((«, v, w), F'(f*, S*, 0, g) [(p, ip, £])y ,y * = ((u, v), ^ )2 (/*, S'*)^ 9)0h,h* 
where 

(0, u), -F 1)2 (/*, = ((u, V, 0), .F'(/*> S*, 0, 0) [<p, V>, 0])y o ,y Q * 



. . (d— S*) 2 , , . i>V> 



+ftmP - + v<p) - k 2 {1 - 3/>y> 



and 



(w,^(f*,g)£) HtH * = ((0,0,w),F'(U,S*,0,g)[0,0,£}) Yo ,Y Q * 



J {w'e + 9 K 2 wi - k 2 (1 - ft)w£} r dr. 



By Theorem 3.1 of |[ABG 99| , when k 2 > 2d 2 the operator T' x 2 > °"* > is bounded away 
from zero (in quadratic form sense.) Hence, if (ip,ip,0 £ ker(^'(/ K , S^, 0, #*)), we take 
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(u, v, w) = (ip, if), 0) and obtain 







> 



((ip, j>, 0), T'(h, S«, 0, g* K ) [<p, V, £])y ,y * 
((<p, Wirtatf*, s *)[<Pi ^])x2,(x2)* 
^(h\\ 2 x + UWx) 



In particular, ip, ip — 0. 

The operator ^(J*,^) = L + g/t 2 where L = — A r — V(r) is a Schrodinger operator 
with potential V(r) = k 2 (1 - f 2 (r)) > and V(r) -> as r — > oo. It is a well-known fact 
in mathematical physics that in dimension two, such operators have at least one negative 
eigenvalue: 

Lemma 4.2 Suppose V : [0, oo) ^ M is continuous, non-negative, V(r) —>■ as r — > oo ; 
and V not identically zero, and define L = —A — V(r) as a self-adjoint operator on the 
space L 2 (M 2 ). Then the ground state energy, 



and is attained at an eigenfunction Uq E H. Moreover, Ao is an isolated, non- degenerate 
eigenvalue, uq G H , and uq > 0. 

The proof follows as an application of the Birman-Schwinger principle in Reed & Simon 
R.S 78|| . We provide an elementary variational proof for the reader's convenience. 




Proof: Let 




Then, 




while 
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(possibly infinite.) Hence, for n = N large but fixed we have J[(u' N ) 2 — V{r)u 2 N ]r dr < 0, 
and hence Ao < 0. Since L is a relatively compact perturbation of —A, Ao is a discrete 
eigenvalue with associated eigenfunction u contained in the form domain of L, H. By 
standard arguments, u > and A is a simple (non-degenerate) eigenvalue. 





By Lemma 4.2 



—g* = inf 

weH-{o} 



j ^K) 2 -(i-/ K 2 ) 



w 



r dr 



J w 2 rdr 



< 0, 



and A = ~^ 2 g* K is the ground state eigenvalue of the Schrodinger operator — A r — k 2 (1 — f 2 ). 
Since Ao is a simple eigenvalue 

dimker(^(/ K ,^)) = dimker(-A r - k 2 (1 - f 2 ) + K 2 g* K ) = 1. 

In conclusion when g = g* the operator J-'{f K , S K , 0, g*) has a simple eigenvalue and the 
eigenvector is of the form (0,0, w K ) with w K the (positive) eigenfunction of T' 3 . 

Finally, we observe that the operator {-§^)T'{f K , S K , 0, g) E L(Yq,Yq), 



r 

((u,v,w),—J r '( y f K ,S K ,0,g)[<f,ij,£})Y ,Y * = J K 2 w£rdr. 



At the eigenvalue g = g* we have 



^-T 1 (L S K , 0, g*) [0, 0, w K ] = k 2 w k & Ran {T\~f Kl S K , 0, g K )). 



Therefore Theorem 1.7 of ||CR 71|| applies, and g* is a bifurcation point for T in Y x M: there 
exists a neighborhood U of (f K , S K , 0, g*) in Y x M, such that the set of non-trivial solutions 
of T(f, S, m, g) = in U is a unique C 1 curve parametrized by ker(.F'(/ K , S K , 0, g*)). 



Remark 4.3 Since T is a smooth (C°°) map, we may calculate various derivatives of the 
bifurcation curve through the normal core solutions at g*. If we parametrize g = 7(i), with 
7(0) = g*, then we follow Crandall & Rabinowitz ||CR 71|| or Ambrosetti & Prodi ||AP 95|| (see 
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Remarks 4.3) to calculate derivatives of j(t) and determine the direction of the bifurcation 
curve locally at g — g*. We obtain that 7'(0) = 0, and 

J f K u*wl + wt. rdr 

7"(0) = ~2 J [J \ " J , 

J wf.r dr 

where u K is obtained from the (unique) solution to the linear system 

^(fn, 3 K , 0, g*) u„ tw*] = - (2K 2 f K w 2 K , 0, 0) 



with u*, w*) -L ker J-'\f K , S K , 0, gr*). By taking the scalar product of the above system 
with (w*, t>*, 0) (and recalling that ^F'(f K ,S K ,0,g*) is positive definite in the complement of 
its kernel) we obtain / f K u*w 2 K rdr < 0, and hence the expression for 7"(0) is indefinite in 
sign. In a joint paper with J. Berlinsky [[ABBG 98|| we present computational evidence that 
solutions bifurcate to the left, to smaller values g < g* K . By standard bifurcation theory (see 
[ |CR 73|| , for example) the direction of bifurcation indicates the stability of the solutions, and 
indeed we observe numerically that the AF core solutions which bifurcate at g* are stable 
(local energy minimizers.) 



4.2 Global bifurcation for (GL)^^ 

We obtain the same abstract bifurcation result for the extreme type-II model, (GL) 
Namely, the value 



inf 

ueH-{o} 



I ( 



; w 



A2 



r dr 



J w 2 rdr 



> 



is a bifurcation point for nontrivial (m > 0) solutions from the (trivial) curve of normal core 
solutions (/oojO, g). But in this case we can make a much more precise statement: 

Proposition 4.4 Let 

E = {(f,m,g) : (f,m) is an admissible solution to (GL)^^ with m > 0}. 

Then C = S U {(/oo; 0, g^)} is a connected C 1 curve, parametrized by g. Moreover for any 
g > , C R {g > go} is compact. 
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As a consequence we have the following exact solvability theorem for (GL)^. 

Theorem 4.5 For g > g*^, the normal core solutions (/oo,0) are the only admissible solu- 
tions of ( GLjoo^g. 

For < g < g^ there is a unique solution with m > 0. This solution is the global minimizer 



The proofs of these two results hinge on the powerful Theorem |3J] and the following com- 
pactness theorem, which will be proven in Section 6: 

Theorem 4.6 Let < a < b. Then the set of all admissible solutions of (GL)^^ with 
g G [a, b] is compact in Z . 

Proof of Theorem |4.4| : Let C be a maximally connected component of C, and suppose 
(fo,m ,g ) G C but (f ,m ,g ) ^ (f^, 0, g*^). Since m = only when g = g*^, we must 
have m > 0. By Theorem [D], (f Q ,m ,g ) is a nondegenerate zero of T in Z x M, so by 
the Implicit Function Theorem there exists a neighborhood U of (/cb m o> <?o) in 2 x J?, an 
interval J = (go — 5, go + 5), and a C 1 function $ : J —>■ Z so that all solutions of JF = in 
U are of the form (<&(g),g) with g G J. 

Let 

<? = sup{g : there exists a solution (/, m, g) G C'} > go- 
Note first that any solution must satisfy 

0< I (m') 2 r dr < I (1 — q)m z r dr. 



and hence g < 1 for any solution with m ^ 0. Since by Proposition [0| C (1 {j > go} = 
C H {(?o < 9 < 1} is compact, there exists a solution at g = g, (/, m, G C. First, we claim 
that m = 0. If not, then by Proposition [2.5| m(r) > for all r > 0, so by Theorem 3.1 , 



(/, m) is a nondegenerate minimum of (GL)§ :00 . By the Implicit Function Theorem argument 
above there exist a C 1 curve of nontrivial solutions through (f,m,g), parametrized by g. 
In particular, we contradict the definition of g is the supremum of all g for solutions in the 
connected component C . Hence rh = 0, as desired. 

Now we show that g = g^. Take a sequence (/„, m n , g n ) G C with g n g, so the above 
arguments imply that f n — /oo — > in X and m n — > in H. Let 



*« = J(l-f n )m 2 n rdr^Q. 
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Then w n = m n /t n solves 

(4.1) - wl - -w' n + g n w n = (1 - f 2 - m 2 n )w n . 

r 

Since 

(4.2) J (K) 2 + gw 2 n ) r dr = J(l - f 2 - n£)v%r dr < J(l - f 2 )w 2 n rdr = 1, 

(by the choice of t n ,) we have ||u>n||# < 1/fl 1 an d we may extract a subsequence (which we 
continue to call w n ) which converges w n — ^ weakly in H and strongly in Lf oc . By the 
strong convergence of /„ — > foo we have /(l — f^w^r dr = 1, so Woo ^ 0, and > 0. 
Passing to the limit in (|4.1|) we have 

(4.3) J (w'^cp' + gw^y - (1 - /^)woo^) r dr = 0, 

for all if G H . This can only occur when g = g^, the ground state eigenvalue of the above 
Schrodinger operator. 

We have just shown that the point (f^, 0, g^) belongs to every connected component of 
C, and hence C is connected. The solution set C is everywhere a C 1 curve: for g > g^ this 
results from the Implicit Function Theorem argument in the first paragraph, and at g^ it 
is a consequence of bifurcation from a simple eigenvalue [|CR 71 1- We now claim that there 



exists exactly one solution in C for every g < g*^. Suppose not, and consider 

D = {g E (0,glo) '■ there exist two distinct solutions (f g> i, fn g ,i), (f g ,2, m g,2) i n C at g.}, 
and g = sup D. 

First, we note that g < g^. To see this we note that the only solution in C with g = g^ 
is the normal core solution, and the bifurcation theorem ensures that the solution set in a 
neighborhood of the bifurcation point (/oo, 0, g£j is a single smooth curve. 

Next, we claim that g ^ D. Indeed, if go G D there exist two distinct solutions 
(/ 90 ,i, iTi go ,i) and (f go ,2, m g ,2) for g = go- By the Implicit Function Theorem argument of the 
first paragraph there exist neighborhoods Ui (of (f go ,i,m g0t i, g )) and U 2 (of (f go ,2, m g ,2, go)) 
in Z x M such that all solutions of T — in Ui, Ui are given by smooth curves parametrized 
by g. In particular, C contains two distinct solutions for g in an interval to the right of go, 
contradicting the definition of go as the supremum. 

Hence g*^ > go G" D, and there exists a sequence —> go for which C contains two 
distinct solutions, (f gh! i, m gkt i), (f gk ,2, m g k ,2)- By Theorem 1 ^. 6|, along some subsequence these 
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solutions converge, and since go ^ D, they both converge to a single solution, (f 9o ,m go ). But 
this contradicts the Implicit Function Theorem argument, which implies that the solution 
set near (f go ,m go , g ) is a single curve parametrized by g. We conclude that the AF core 
solutions are unique for each g £ (0, p^). 





4.3 Behavior for g — > 0, k < oo 



For the problem (GL) K)S we do not have the strong information provided by Theorem |3J 
which determines the global structure of the solution set, and hence we cannot make the 
same elegant conclusion about the uniqueness of AF core solutions. However we may still 
say something about the global structure of the continuum bifurcating from the normal cores 
at g — g*. When n 2 > 2d 2 we may apply the Global Bifurcation Theorem of Rabinowitz 



Ra 71|| to conclude that the continuum S K of zeros of J-(f, S,m } g) = with m > is 
unbounded in the space Y x M. (Note that E K cannot contain any other eigenvalues of the 
linearization about the normal core solutions, as is easily seen from the calculations (|4.1|) - 
( |4.3|) above.) In the next section we will prove the following a priori estimate, which has as 
a direct consequence the fact that S K can only become unbounded as g — > + : 



Theorem 4.7 Let d,n be fixed. For any compact interval J G (0, oo) there exists Co = 
Co(K,d,J) > such that every admissible solution (f,S,m) of (GL) K ^ g with g E J satisfies 
\\(f,S,m)\\ Y <C. 

Let us now concentrate on this loss of compactness in the continuum S K as g — > + . We 
prove: 

Theorem 4.8 For any sequence of (absolute) minimizers (f g ,S g ,m g ) £ Y with g —>■ + we 
have f g —>-0in Xi oc , S g — > locally uniformly, and m g — > 1 in Hi oc . 

Fix k £ M, and for any g > consider a minimizer (f g ,S g , m g ) £ Y of £ K ^ g . 

Lemma 4.9 

£ K , g (f g ,S g ,m g ) — >0, as g^O. 

Proof: We will show that for any e > there exist g e > and H radial functions 

(f e , S £ , m £ ) £ Y such that < £ K>g {f e , S £ , m £ ) < e, for any g < g £ . 
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For a fixed p > 0, we define 



if r < p, 

VI' ) ~ \ ln(l/p) 



we consider 



and 



if r > p 



7T 7T 

f p (r) = cos(u p (r)-), m p (r) = sm(u p {r)-), 



S P (r) 



0, ifre(0,§) 
d, if r G (p, oo 



2/' 



A direct computation shows that 



_ , _ > C^j IT hC a 

S K!9 (f p , S p , m p )<- + ^— + —gp 



for any g > 

p2 41np e ^ 2 



For a given £ > 0, we choose a p e such that ^ + jf 1 — < |, and a g e = g £ (p £ ) for which 



< I' Le - 9e < |^4- Th en, S g ,m g ) < £ K ,g(f Ps ,S Pe ,m Ps ) < e, for any g < g £ . 





Proof of Theorem [4.8| By Lemma |4.9| each term in the energy tends to zero as g — > 0. 
First, note that j(S' g /r) 2 r dr — ► combined with (1.4) in | BC 89|| implies that 



(4.4) S g (r)/r — > uniformly. 

For any Rq > 0, we then have 



* (rf-^) 2 , 2 , r R ° (d-S q ) 2 e2 , 
o(l) = ^-f g rdr> K - ^f g rdr 



r* Jo r 



r R d 2 
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In particular, f 9 — > in Lf oc , X\ oc . Finally, by the reverse triangle inequality, 



V2 l 



2 ([o,floD 



> 



K 

71 



1 - "^^([O^o]) - ll/ 9 2 |U 2 ([o,fl ]) 



> ^ll 1 -^ 9 IU 2 ([0,fl ]) +0(1), 



where we have also used < f g < 1, < m g < 1, and f g — > in Z^ oc . In conclusion m g — > 1 
in L; 2 oc and in fact in -ff/ oc , since J(m' g ) 2 rdr —>■ by the energy estimate. 



5 The limit ^ — > oo 

In this section we show that the problem (GL)^^ arises as a limiting case of (GL) KjS as 
k — > oo. For any solution (/ K , S K , m K ) of (GL) Kjff , define 

(5.1) f K ( r ) = f K (£j, S n (r) = S K (~), Mr)=m K (~y 

We prove: 

Theorem 5.1 Lei S K , m K ) be any family of solutions of ( GL) Kj9 for k > 0, and (/ K , S K , rh K ) 
defined as in ( \5.1\ ). For any sequence n n — > oo ; there exists a subsequence and a solution 
(/oo, moo) of ( GL)^ so that (as n nk -> oo J f Kn - f^ in X , rh Kn - -> m if, and 
S Kn — > locally uniformly. Moreover: 

»■ #0 > fi^, ^en m K -> 0; 

M. 7/m K ^ /or a// /arge K and g ^ g*^, then lim^oo m K = > 0. 
As a simple consequence of the uniform convergence of f K — > /^ we have the following 
Corollary 5.2 
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Remark 5.3 This implies that the bifurcation diagram for (GL) Kj9 with k very large should 
strongly resemble the very precise image given for (GL)^^ by Theorem |4.5| . In particular, 
for any fixed g > (GL) K S cannot have solutions S Kj9 , m^ g ) with m Kg > for k 

large. 

Simple calculations using the energy £ Kg show that infy£ Kg ~ Inn, and hence we re- 
quire require energy-independent estimates for our solutions (f K , S K ,m K ). To obtain these 
estimates we begin with a simple version of the celebrated Pohozaev identity. This identity 
will also be essential for proving the a priori estimates used in the bifurcation analysis in 
the previous section. 

Proposition 5.4 For any finite energy solution (f,S,m) of (GL) K ^ g we have 

m 2 r dr + — J (1 — f 2 — m 2 ) 2 r dr = J — r dr. 

For any finite energy solution (f,m) of (GL)^^ we have 

f If d 2 

g J m 2 r dr + — J (1 — f 2 — m 2 ) 2 r dr = —. 

Proof: We multiply the first equation in (GL) Kj9 by f'(r)r and integrate r dr to obtain: 



K 

~2 



(1- f 2 -m 2 )(f 2 )'r 2 dr 



(d-S)S'f 2 dr 



— ) S' r dr 

r 



J ^ S 'r)'dr = J (y) rdr, 



using the equation for S(r), and integrating by parts whenever necessary. We also multiply 
the third equation in (GL) KS by m'(r)r and integrate rdr to obtain: 

— J {1 — f 2 — m 2 ){m 2 )'r 2 dr = ^— J {m 2 )'r 2 dr = —gn 2 J m 2 r dr. 
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Together, 



f ( S'\ 2 r n 2 r 

J ( — j r dr = gn 2 J m 2 r dr + — J (1 — f 2 — m 2 ){m 2 + f 2 )' r 2 dr 

= gn 2 J m 2 r dr + — J (1 — f 2 — m 2 ) 2 r dr. 
For the case k = oo we proceed in the same way, except the equation for / yields 



l_ / 2_ m 2 ) (^] r 2 dr 



2 ' 



The calculation then continues as above. 



Proof: Step 1: Bounding the sequence. 



From the Pohozaev identity (Proposition |5.4j ) and Lemma 4.2 of | |BC 89| | after rescaling we 
have: 



d 2 > 



— rdr 
r I 



(5.2) 



i I rdr 



(5.3) 



g j m 2 K rdr + ^ / (1 - f 2 - m 2 K ) 2 r dr. 



Using (|5.2|) and Lemma 1.2 (ii) in PC 89|| we have 



(5.4) 



sup 

r€[0,oo) 



and hence S K — » locally uniformly. From (|5.3| ) we obtain the uniform bound ||m K || 2 < C 
(depending on g, which we assume is fixed.) From the equation for m K , after a change of 
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scale, we obtain: 



■m K ) + gm K 



rdr — / (1 — f 2 — rh 2 .) r dr < / m 2 K r dr < C 



and therefore ||m K ||i/ < C uniformly in k. 

Recalling Proposition any solution satisfies < rh K (r) < 1, and we may conclude 
that ||m rc || g < ||rri K || 2 < C for all p e [2, oo]. By the triangle inequality, 

||1-/ K 2 || 2 < ||l-A 2 -m2|| 2 + ||^|| 2 



and hence we obtain 



< \d\ + C, 
(l-f K ) 2 rdr< [{l-flfrdr<C 



(since f K > 0.) 

Choose a function rj G C°°(M) with 



rj(r) 



1, ifr<2, 
0, if r > 3, 



and < rj(r) < 1 for all r. Using r] 2 f K as a test function in the weak form of the rescaled 
equation for f K , 



2 d2 "2 



r dr 



(i-f 2 -Of 2 v 2 -vv'fJL 



r dr 



< 



^(l-/ K 2 ) 2 + ^ 4 + ^ 2 (/«) 2 + 2/^') 2 ' 



r dr 



1 



< C + ~J V 2 (f K ) 2 rdr. 
Absorbing the last term back to the left hand side, 



(5.5) 



ULY + p 2 



rdr < C. 
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Now choose another smooth function fo with 



/o(r) = { I 



0, if r > 1, 
if r > 2, 



and < / (r) < 1. Note that with this choice fo+r] 2 > 1. We use (f K — 1)/q as a test 
function in the equation for f K to obtain 



(£) 2 + ^(/ 2 -i) 2 



(1 - # - mi)Ufn - l)/ 2 - (/« - l)/«/o/o 



r dr 



4(^- i )/o 2 + 



S K (2d S K 



2 (/« l) 2 /o 



r <ir 



< tfjf (/„-l)|£|/ rdr 



+ 



d 2 4d - 

_(!-/„) + _(/„-!) 



r dr 



< \ liftffirtlr + CWf,. 



+ 



§ + (A 2 - 1) 4 



r dr 



(Note that in the first line, the first integrand is non-positive.) In conclusion, 



(5.6) 



d 2 



r dr < C. 



Now define u K = f K — f e X. Then from (|5.5|), (|5.6|) we obtain: 



d 2 ' 

K) 2 + u 2 K + 



r dr < 2 



2 2d 2 - 
rdr+l —z-K r dr 
o r 2 
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+ J 2 [d 2 + l){f K -lf 



rdr + C 



< c. 

In other words, u K is uniformly bounded in X, and we may extract weakly convergent 
subsequences u n = u Kn — 1 u* (in X), m n = m Kn — > m* (in H). 

Step 2: Strong convergence. 

We next show that the sequences u n , m n converge in norm. Let f n = fo + u n and S n = S Kn . 
First note that 



(1 - fn ~ m 2 n )m n - (1 - /p - ^)m p = (1 - f 2 )w - (ml + m n m p + m 2 p )w + (f 2 - f 2 ) 
Hence, using compact embeddings of X, H into L q for 2 < q < oo, 



m. 



(K-mp)) +g(m n -m p ) 



r dr 



= J [(1 - /*) - (m* + m n m p + m 2 p )(m n - m p ) 
+ Jifn + fp)(u n - u p )m p (m n - m p )rdr 



(m n — m p ) 2 r dr 



= o(l). 

Therefore, m n — > in norm. 

We proceed in the same way with u n : 



/{k-o 2 + 



„.2 „.2 •'P 



(5.7) 

Now we expand, 



/ {(!-/«- - (1 - fp ~ m 2 p )f p ] (u n - u p )rdr 



(d-S n ) 2 (d - S p ) 2 

9 In 9 Jp 



" -^(2d-5 B ) 



(« n - M p ) 
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^(2d-S p )-^(2d-S n 



Now we take each term separately: 



/ -%(2d- S n )(u n -u p ) 2 rdr < sup \S n \ / 2d 
Jo r z r6[0 ,i] Jo 



Un Up) 7 



since S n — > locally uniformly, and U n £1X6 uniformly bounded. 



poo 

-^r(2d - S n )(u n -u p ) 2 rdr < sup 



5n 




r 





2d(u n — u p ) 2 rdr^O, 



by (|5.4j). Choose r > so that / r °°r 2 dr < e 2 /4, and k sufficiently large so that 



drp sup 



5*,, 



\Un - WplU < 2' 



Then, 



fro S S If 

I f p -j(2d- S p )(u n - Up) r dr < dr 2 sup —\ 

JO r Z 0<r<r r y Jo 



ro { u - up) 2 

— - — r dr < e/2, 



and 



/ fp^{2d-Sp){u n -u p )rdr <2d 2 

Jr r z 



r°° dr 


1/2 


[/ 


Jro r 3 _ 




Jo 



(u n — Up) 2 r dr 



1/2 



< e/2. 



We return to (|5.7|) , and substitute the above estimates: 



{u' n - u'p) 2 + ^(u n - Up) 2 \ r dr + <>( I ) 
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= J {(1 - fn ~ m l)fn ~ (1 - fp - m 2 p )f p } (u n - u p )rdr 
= J {(1 - 3/ 2 )K - Up) 2 - 3f (u 2 n - u 2 p ){u n - Up) 

-{ul - u 3 p )(u n - Up) - f (m 2 n - m 2 p )(u n - u p ) 
-m 2 n (u n - Up) 2 - u p (m 2 n - m 2 p )(u n - u p )} rdr 

= -2 J \u n - u p ) 2 rdr + o(l), 

where we use the facts that m n — > strongly in H, u n is bounded in X, and u n — > in 
Lf oc . In conclusion, the subsequence u n — > strongly in X. 

Step 3: Determining when = 0. 

Since all solutions of (GL)oo i9 with g > have = 0, we have m K — > when g > g^. 
On the other hand, suppose m K > for all sufficiently large k, but m K — > 0. By uniqueness 
of the normal core solution, f K — > /qq, the unique solution of 



Ajoo + = (1 - £>)/o 



Let 



t« = J(l-f*)rh 2 K rdr^0, 



and set u> K = rh K /t K . Then 



-A r w K + ^ K = (1 - f K - m 2 K )w K . 



Since 



J [«) 2 + gw 2 ] rdr = J(l-f 2 - m 2 K )w 2 K r dr < 1 



(by the choice of t K ,) the bound < l/g results. We extract a subsequence (which we 

still denote by w K ) with w K — 1 Woo weakly in H. Note that u>oo > 0. By the choice of t K , the 
uniform convergence — > /oo) and the Lf oc convergence of w K — > we have: 



/ [(i - - + (A 2 - LX + (i - 



r dr 



= l + o(l). 
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In particular Woo ^ 0. By weak convergence we may pass to the limit in the equation for 
w K , and hence to^ is a nontrivial non-negative solution of 

-A r Woo + QWoo = (1 - /^)Woo- 

This can only occur when g = g^. 



This completes the proof of Theorem |0 

<> 



6 Estimates and existence 

In this section we derive the technical estimates which were needed in our analysis of the 
bifurcation problem in Section 4. We also provide the details of the proof of existence of 
minimizers of the energies £ Kt9 and £oo, g - 



6.1 A priori estimates 

We may now prove a priori estimates for the solutions of our system (GL) KS , Theorem |4.7| , 
as well as the compactness result for solutions of (GL)oo j9 (both theorems as stated in the 
previous section.) Note that both theorems are stated for all solutions, not only energy 
minimizers, and hence we will use our Pohozaev identity (Proposition |5.4j) to obtain energy 
independent estimates. As before, we denote by f K , S K a normal core solution at k, and 
u = f - f K , v = (S - S K )/r. 



By the Pohozaev identity and Lemma 4.2 of [BC89] we have 



(6.1) 



K 



gm 2 + -(1 - f - 'in 



2\2 



r dr 



'S'Y J d 2 

— r dr < — . 
r ~ 2 



In particular, we obtain 

J {I- ffrdr <C + C/g, Jm 2 rdr<C/g 
with constant C depending on k, d. From the first estimate we obtain 

\\uh < Win + u - 1|| 2 + || A - 1|| 2 < C + C/g. 
The equation for m together with the second estimate gives: 

/ [m') 2 + K 2 gm 2 r dr = k 2 / (1 — f 2 — m 2 )m 2 r dr < n 2 I m 2 r dr < C/g. 
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In particular, ||j7i||h < C, IMh < C, and the constant depending on k, d may be chosen 
uniformly for g G J. 



Using the right half of (|6.1| ) we have 



(6.2) 



d 2 f ( S'\ 
— > / I — r dr 
2 - J r 



r dr. 



Since 



2 / -JLi — '-rdr 



< 2 



r 



+ 



[rv 



r dr, 



and 



[rv 



r dr 



(v'Y + 



r dr. 



we may conclude from (|6.2j) that 



(6.3) 



.A 2 



rdr < C, 



with constant depending only on d. From the embedding properties of X, Lemma |2.1| , we 
conclude that ||f ||oo < C. 

We now use v as a test function in the weak form of the equation for S to obtain an 
estimate: 



(6.4) 



d-S 



f 2 v r dr 



S' 

— (rv) 1 r dr 



< 



's'Y . i 

— rdr + - 



(v'Y + 



rdr < C. 
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On the other hand, expanding the left-hand side of flB.4|), 



(6.5) 



'd-S" 



f 2 v r dr 



d — S K 



(f K + u) 2 vr dr 



d — S, 



[f K + 2f K u + u }vr dr — \ 2f K uv r dr 



v 2 fir dr — I v l u l r dr. 



.22 



To bound the term / v f 2 rdr, we need to evaluate the other terms: 



j K vr dr 



< 2 



d — S K 



1 f 1 

fl r dr + - J f 2 v 2 rdr <C + - / f 2 v 2 r dr, 



d — S K ~ 
j K uv r dr 



- n 2 



<-I 


d- S K 






r 



f*rdr+ \\v\UuU < C 



d S K o 



u vr dr 



If. 1 1 9 1 1 no 1 

<2«ll* + 2 



d- S K 



u rdr, 



f K uv 2 r dr 



<8|H|iWL + ^ (f~yrdr<C+l I f 2 K v 2 rdr, 



v 2 u 2 rdr < \\v III^IImIIo < C. 



2 ll„.l|2 



Hence, 
(6.6) 



-Jv 2 f 2 rdr<C+ 1 - 



d — S K 



u rdr 



Finally, we use u as a test function in the weak form of the equation for /. Recalling the 
definition of f K as a normal core solution, we expand and cancel terms to arrive at: 



(6.7) 



,'\2 



' d - S K 



r dr 



2 ( <LJk\ ( ; k + u)uv _ ( / k + u)uv i 



r dr 
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+k 2 J [(1 - 3f 2 )u 2 - 3f K u 3 -u 4 - m 2 fu 
Each term on the right hand side may be controlled as follows: 



r dr. 



m 2 fur dr 



1 f 1 
< — J m 4 + - / u 2 r dr < (' . 



(1 - 3f 2 )u 2 rdr 



<3\\u\\l<C, 



'd-S, 



— ] f K uv r dr 



<\\u\\l\\v\\l+ I flrdr<C, 



f K u 3 r dr 



d S K \ 2 



u vrdr 



<6|Hll|H|L + i/(^) u 2 rdr 



f K uv 2 r dr 



< \\Hl\Hl + ~ / fyrdr <C+U (i-A) u 2 rdr. 



where in the last estimate we apply (|6.6|) . Using (|6.?1 ) we have 



r dr < C. 



Consequently, < C. Returning to (|6.6|) , it follows that 



/ 



fyrdr<C+~J (*—?*\ u 2 rdr<C, 



and hence (|6.3| ) yields \\v\\x < C- This concludes the proof of Theorem [O 
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An analogous result may be proven for solutions of (GL)oo i9 : 

Theorem 6.1 Let d be fixed. For any compact interval J G (0, oo) there exists Cq = 
Co(d, J) > such that every admissible solution (f,m) of (GL)^^ with g G J satisfies 
\\(f,m)\\ Zo <C. 

The proof of Theorem |6.1| is similar to (and simpler than) the previous one, and is left to 
the reader. 



6.2 Compactness 



Here we prove Theorem [4.6| , which asserts that the family of solutions to (GL)oo iS with g 
bounded away from zero is a compact set. The same result holds for (GL) Kj9 , although the 
proof is more complicated due to the additional terms involving S(r). 

Proof of Theorem |4.6j : Suppose f n = foo + u n , m n , g n are a sequence of solutions of 
(GL)oo,g n with g n G [a, b}. By the Theorem we have ||u n ||x, ||^n||j? < C, and hence we may 
extract a subsequence with u n — ^ u, m n — ^ rh, and g n — > g G [a, b\. Then we have 



(G.i 



d 2 

K - U 'kY + ^(Ur 



Uk 



r dr 



I 1 ~ fn - m l)fn ~ (1 - fk - m l)fk K - u k ) r dr 



(6.9) 



(m' n - m' k ) + g(m n - m k y 



r dr 



(! - fn~ m l) m n ~ (1 - fk ~ m l) m k\ K - m k )rdr + o(l). 

We now expand the two right-hand side terms. First, we use the embedding properties of 
X, H and the fact that < f n < 1 for any solution to show: 



;i - f n - ml)f n - (1 - fl - m 2 k )h 



{U n - Uk) 



Vl{u n - u k ) 2 r dr + / (1 - f^)(u n - u k ) 2 r dr 



-2 / [/oo(m„ + u k ){u n - u k ) 2 - f n (u n + u k ){u n - u k ) 2 ]r dr 
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- / [fn(m n + m k ){m n - m k ){u n - u k ) + (u k + m k )(u n - u k ) }rdr 



2~f 00 {u n -u k ) 2 rdr + o{l). 



Applying the above estimate to (|6.8|) we have 



K - u k f +{-,+ 2f~A K - u k ) 



r dr — > 



as n, k — ■> oo, so w n — > w in norm on the space X. 
Similarly, we estimate 

(1 - fn ~ m l) m n ~ (1 - fl ~ m 2 k )m k (m n - m k ) r dr 



y (! - fL){ m n ~ m k ) 2 rdr - J ftf^u^rrin - m k f 



+2f 00 m k (m_m k )(u n - u k ) + u 2 n (m n - m k ) 2 ]rdr 



- / [m k (u n - u k ) {m n -m k ) + (m n - m k ) (m n -m k )]rdr 



= o(l). 

Therefore, ( |6.9|) implies that m n —>■ fh in H . By passing to the limit in the weak formulation 
of (GL) 5ni00 we easily obtain that (f,m) solve (GL)^ )00 , and hence the specified solution set 
is compact. 



6.3 Existence 

Let (u n , v n , m n ) be a minimizing sequence for I Kj9 , so (/„, S n , m n ) = (f + u n , S + rv n , m n ) 
is a minimizing sequence for To prove Theorem [T2] we first observe that the energy 

S Kt9 is a sum of positive terms, and hence each is individually bounded. In particular, m n is 
uniformly bounded in H. 
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Now we must estimate u n . First, note that £ K ,g{\fn\i S n , m n ) = ^K, g (fn, S n m n ), and so we 
may assume that our minimizing sequence satisfies f n { r ) > for all r. Next, we observe 

(6.10) ||1 - /„ 2 || 2 < ||1 - f 2 n - m 2 || 2 + ||m 2 || 2 < ||1 - f n - m\ || 2 + C. 

Hence we conclude that 



> 



r 



+ 



ft 



• T si-./;) 2 ^ -//•. 



The right-hand side of the above inequality is the free energy of conventional Ginzburg- 
Landau vortices studied in ||ABG 99| . The boundedness of ||w n |U, Iknllx then follows from 



the argument of Proposition 4.2 of |[ABG 99| . We may then pass to the limit in £ K g via 
lower semicontinuity of the norms and Fatou's Lemma. 



To prove Theorem |2]^ let (u n , m n ) be a minimizing sequence for / m in I x H, so 
(fn, m n) = (Zoo + u n , m n ) is a minimizing sequence for £oo, 9 - Choose r g > 1 so that ^ < |. 
Then 



£oo,g(fm Win) — 



d 2 ~ ' 

I I \2 i 2 u x2 



r dr 



I [{<? + 0/ - ^VO™ 2 + -Ml + < ~ 1) 



+ ^(/ 2 + ™ 2 -l) 2 + ^(l-/^)kdr 



> 



<- / \2 , 9 2 



r (ir 



o r 



d 2 ~ 

— /£, r rfr + 



r/ 2 ~ r/ 4 ' 

— ri - f 2 ) - — 

r 2\ •'oo' 1 2r 4 



r <ir, 



where we have used the elementary bound ax + x 2 /2 > — a 2 /2. In particular, £00,3 is 
bounded below and the minimizing sequence has ||m n ||H < C uniformly in n. By the 
Sobolev embedding, we also conclude that ||?«n||p < C p for all p G [2, oo). 
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Now we must estimate u n . As above we note that Soo^d f n \,m n ) = £oo,g(fn, m n) , and so 
we may assume that our minimizing sequence satisfies f n { r ) > for all r, and the bound 
( 16.101) holds. Note that we also have: 



(6.11; 



K|| 2 < ll/oo - 1|| 2 + 111 - fnh <C+\\l-f, 



n 2- 



By the estimate on m n , ( |6.10| ), and (|6.11| ) we now have 



C > 



d 



1 



2\2 



(fnY + ^ifn-fD + ^-fn) 



r dr 



r dr 



U' n ? + (f'J 2 + -A + \( l - fnf + 2(1 - fDlooUn 



r dr 



> 



d 2 1 
{ / \ 2 1 2 1 2 
K u„) + — u„ + -u n - 



1 9 



r 



2 « 



4 



1 /oo) O foo u n 



r dr — C 



> 



r/ 2 1 

/ / \ 2 2 2 



r 



2 n 



r dr — C. 



In conclusion ||w n ||x < C- We extract a subsequence for which both u n — ^ -u an d ^ n - "io 
weakly in X, H respectively, and pointwise almost everywhere. 

By semicontinuity of the norm, Fatou's Lemma (for the positive terms) and the l? rloc 
convergence of u n — > u we can pass to the limit in 



Ioo,g{uo,m ) < liminf Ioo(u n ,m n ) = inf 

' a n-»oo XxH 



So the infimum of 1^ is attained. 





S. Alama, L. Bronsard, T. Giorgi 



42 



References 

[ABG 99] S. Alama, L. Bronsard and T. Giorgi, Uniqueness of symmetric vortex so- 
lutions in the Ginzburg-Landau model of superconductivity, to appear in 
Journal of Functional Analysis. 

[ABBG 98] S. Alama, A. J. Berlinsky, L. Bronsard and T. Giorgi, Vortices with antifer- 
romagnetic cores in the S0(5) theory of high temperature superconductivity, 
preprint |cond- mat / 9 8 1 2 2 83| . 

[ABKZ 97] D. P. Arovas, A. J. Berlinsky, C. Kallin and S.-C. Zhang, Superconducting 
Vortex with Antiferromagnetic Core, Phys. Rev. Lett., vol. 79 (1997), pp 
2871-2874. 

[AP 95] A. Ambrosetti and G. Prodi, "A Primer of Nonlinear Analysis" , Cambridge 

Studies in Advanced Mathematics vol. 34, Cambridge University Press, Cam- 
bridge, 1995. 

[BC 89] M. S. Berger and Y. Y. Chen, Symmetric vortices for the Ginzburg-Landau 

equations of superconductivity and the nonlinear desingularization phe- 
nomenon, Jour. Functional Analysis, vol. 82 (1989), pp 259-295. 

[BMR 94] H. Brezis, F. Merle, and T. Riviere, Quantization effect for —Au = u(l — \u\ 2 ) 
in R, Arch. Rat. Mech. Anal., vol. 126 (1994), pp35-58. 

[CEQ 94] X. Chen, C. M. Elliot, and T. Qi, Shooting method for vortex solutions of 
a complex-valued Ginzburg-Landau equation, Proc. Royal Soc. Edinburgh, 
vol. 124A (1994), pp. 1075-1088. 

[CR 71] M. Crandall and P. Rabinowitz, Bifurcation from simple eigenvalues, J. Func- 

tional Analysis vol. 8 (1971), pp. 321-340. 

[CR 73] M. Crandall and P. Rabinowitz, Bifurcation, perturbation of simple eigenval- 

ues and linearized stability, Arch. Rational Mech. Anal. vol. 52 (1973), pp. 
161-180. 

[GT 83] D. Gilbarg and N. Trudinger, "Elliptic Partial Differential Equations of Sec- 

ond Order", Springer- Verlag, Berlin, 1983. 

[JT 80] A. Jaffe and C. Taubes, Monopoles and Vortices, Birkhauser, Boston, 1980. 



S. Alama, L. Bronsard, T. Giorgi 



43 



[Gu 98] S. Gustafson, "The Stability of Abelian Higgs Vortices," talk presented in the 

session MS75 (Nonlinear Partial Differential Equations) at the SIAM Annual 
Meeting, Toronto ON, Julyl3-17, f 998. 

[LL 94] E. Lieb and M. Loss, Symmetry of the Ginzburg-Landau minimizer in a disc, 

Math. Research Lett., vol. f (1994), 701-715. 

[M 95] P. Mironescu, On the stability of radial solutions of the Ginzburg-Landau 

equation, Jour, of Funct. Analysis, vol. 130 (1995), 334-344. 

[OS 97] Yu. N. Ovchinnikov and I. M. Sigal, Ginzburg-Landau Equation I. Static 

Vortices, "Partial differential equations and their applications (Toronto, ON, 
1995)", pp. 199-220, CRM Proc. Lecture Notes, 12, Amer. Math. Soc, Prov- 
idence, RI, 1997. 

[P 80] B. Plohr, The Existence, Regularity, and Behavior of Isotropic Solutions of 

Classical Gauge Field Theories, Department of Physics, Princeton University, 
PhD dissertation, 1980. 

[P 81] B. Plohr, The behavior at infinity of isotropic vortices and monopoles, J. 

Math. Phys., vol. 22, (1981), 2184-2190. 

[Ra 71] P. Rabinowitz, Some global results for nonlinear eigenvalue problems, J. 

Functional Analysis, vol. 7 (1971), pp. 487-513. 

[RS 78] M. Reed and B. Simon, "Methods of Modern Mathematical Physics. IV. 

Analysis of operators," Academic Press, New York-London, 1978. 

[S 94] 1. Shafrir, Remarks on solutions of — Au = (1 — |w| 2 )w) in M 2 , C. R. Acad. 

Sci. Paris, tome 318, Serie I (1994), 327-331. 

[St 77] W. Strauss, Existence of solitary waves in higer dimensions, Comm. Math. 

Phys, vol. 55, (1977), 149-162. 

[Str 90] M. Struwe, "Variational Methods," Springer, Berlin, 1990. 

[T 96] M. Tinkham, "Introduction to Superconductivity," 2nd edition, Mc Graw- 

Hill, New York, 1996. 

[Z 97] S.-C. Zhang, A unified theory based on SO(5) Symmetry of superconductivity 

and antiferromagnetism, Science, Vol. 275 (1997), pp 1089-1096 



